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Abstract 

In recent investigations, it has been found that conservation laws generally lead to precision limits 
on quantum computing. Lower bounds of the error probability have been obtained for various logic 
operations from the commutation relation between the noise operator and the conserved quantity 
or from the recently developed universal uncertainty principle for the noise-disturbance trade- 
off in general measurements. However, the problem of obtaining the precision limit to realizing 
the quantum NOT gate has eluded a solution from these approaches. Here, we develop a new 
method for this problem based on analyzing the trace distance between the output state from 
the realization under consideration and the one from the ideal gate. Using the mathematical 
apparatus of orthogonal polynomials, we obtain a general lower bound on the error probability 
for the realization of the quantum NOT gate in terms of the number of qubits in the control 
system under the conservation of the total angular momentum of the computational qubit plus the 
the control system along the direction used to encode the computational basis. The lower bound 
turns out to be more stringent than one might expect from previous results. The new method is 
expected to lead to more accurate estimates for physical realizations of various types of quantum 
computations under conservation laws, and to contribute to related problems such as the accuracy 
of programmable quantum processors. 
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I. INTRODUCTION 



Recently, there have been extensive research efforts to explore whether fundamental phys- 



ical laws put any constraints on realizing scalable quantum computing. Soon a: 



ter the dis- 
^ that the 



covery of Shor's algorithm it was pointed out by several physicists 
decoherence, the exponential decay of coherence in time, caused by the coupling between 
a quantum computer and the environment would cancel out the computational advantage 
of quantum computers. To overcome this difficulty, quantum error-correction was proposed 
[s], [g] , and the subsequent development has established the so-called threshold theorem: if 
the error caused by the decoherence in individual quantum gates is below a certain constant 
threshold, it is possible in principle to efficiently perform an arbitrary scale of fault-tolerant 
quantum computation with error-correction [3]. Thus, the error-correction reduces, in prin- 
ciple, the scalability problem to the accuracy problem requiring individual quantum logic 
gates to clear the error threshold, though being still quite demanding. 

In general, decoherence in quantum computer components can be classified into two 
classes: (i) static decoherence, arising from the interaction between computational qubits, 
typically in the memory, and the environment, and (ii) dynamical decoherence, arising from 
the interaction between computational qubits, typically in the register, and the control 
system of gate operations 8]]. The static decoherence may be overcome by developing 
materials with long decoherence time. On the other hand, the dynamical decoherence poses 
a dilemma between controllability and decoherence; the control needs coupling, whereas 
the coupling causes decoherence. Thus, even if the interaction with the environment is 
completely suppressed, the error caused by the dynamical decoherence still remains. Clearly, 
if the control system is described classically, there is no decoherence. However, this never 
happens in reality with finite resources. 

Barnes, Warren 9|, Gea-Banacloche lOj, van Enk, and Kimble 11| have been focused 
on the atom-field interaction between atom qubits and control electromagnetic fields, and 
shown that, when the control field is in a coherent state, the gate error scales as the inverse 
of the average photon number. In contrast to those mo del- dependent approaches, one of 



the authors 



12l | explored the physical constraint on the error caused by dynamical decoher- 



ence generally imposed by conservation laws and obtained accuracy limits by quantitatively 
generalizing the so-called Winger- Araki-Yanase theorem isl, [l^ : observables which do not 
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commute with bounded additive conserved quantities have no precise and non-disturbing 
measurements. It is natural to assume that conservation laws are satisfied by the interaction 
between the qubit and the external control system. If the control system were to be com- 
pletely described as a classical system, the conservation law would not cause any confiict 
in realizing a unitary operation on the computational qubit, since the classical interaction 
causes no decoherence and yet conserves the (infinite) total quantum number. However, 
in reality, the interaction may cause decoherence and the time evolution operator on the 
composite system is limited to one commuting with the conserved quantity. Under these 
conditions, the accuracy of the realized gate operation generally depends on the kind of 
gate being considered. It has been shown that the SWAP gate can be realized in principle 
without error However, the controlled NOT gate and the Hadamard gate have lower 
bounds of the error probability that scales as the inverse of the size of the control system, 
as follows. 

The impossibility of precise and non-disturbing measurements under conservation laws 
was generalized to an inequality for the lower bound of the sum of the noise and the dis- 



turbance of measuring process under a conservation law 



This inequality leads to a 



general lower bound for the error probability of any realization of the controlled-NOT gate 



under conservation laws 
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151 ]. For single-spin qubits controlled by the A^-qubit control 



system, the angular momentum conservation law leads to the minimum error probability 
(4A^^)~^ |l^- Thus, assuming the threshold error probability 10^^ - 10^^ , a two- qubit 
unitary operator needs to be realized by an interaction with more than 100 qubit systems, 
suggesting the usefulness of schemes based on multiple-s pin encoded qubits such as the uni- 
versal encoding based on decoherence-free subspaces 



15 



16 



17j. In bosonic controls, such as 
electromagnetic fields in coherent states, the minimum error probability amounts to {lQn)~^ 



12l |. where n is the average number of photons. The above result also leads to a conclusion 



that in any universal set of elementally logic operations there is at least one logic operation 
that obeys the error limit with the same scaling as above 121 ]. 

On the other hand, without assuming the non-disturbing condition the lower bound for 
the noise in arbitrary measurements under arbitrary conservation laws was derived from 
the commutation relation for noise operator and the conserved quantity [l^ or simply from 



the universal uncertainty principle 



19| : see Refs. 



21 



221 ] for the universal uncertainty 



principle. This inequality also leads to a general lower bound for the error probability of the 
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realization of the Hadamard gate that amounts to the minimum error probabihty (4A^^) ^ for 



any A^-qubit control system and (16n) ^ 
state with average number of photons n 



br any electromagnetic control field in a coherent 



19| . Gea-Banacloche and one of the authors 23 1 



compared the above result for electromagnetic control fields with the previous result obtained 
by Gea-Banacloche 10|] for the Jaynes-Cummings interaction, and it was concluded that the 
constraint based on the angular momentum conservation law represents an ultimate limit 
closely related to the fluctuations in the quantum field pha se. The use of the Jaynes- 
Cummings model in the above model-dependent approach 
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llj was questioned by Itano 

in 



26| 



24| and subsequently Silberfarb and Deutsch j25|] justified the Jaynes-Cummings mode 
the limit of small entanglement; see also replies to Itano by van Enk and H. J. Kimble 
and by Gea-Banacloche [27]. The above consistency result between the model-dependent 
and model-independent approaches enforces the validity of the use of the Jaynes-Cummings 
model and substantially clarifies the whole situation. 

The above methods for deriving conservation-law-induced quantum limits for quantum 
logic operations are also applicable to the Toffoli gate and the Fredkin gate to obtain similar 
lower bounds. However, the problem of obtaining the precision limit to realizing the quantum 
NOT gate has eluded a solution from these approaches, and hence the problem has been 
open as to how the minimum error for that gate scales with the size of the control system. In 
this paper, in order to solve this problem we devise a new method of deriving the precision 
limit, and show that there exists a non-zero lower bound, which indeed scales as the inverse 
size of the control system, of the error probability for the quantum NOT gate. 

Our formulation has various common features with the formulation of programmable 
quantum processors 



28l . |29| . |30|, in which a set of unitary operators is to be realized by 
selecting a unitary operator on the composite system, the system plus the ancilla, and 
by selecting a set of ancilla states, whereas in our problem a single unitary operator is 
to be realized under a conservation law by selecting a unitary operator on the composite 
system satisfying the conservation laws and by selecting a single ancilla state. In previous 
investigations the accuracy of programmable quantum computing has been measured by the 
so-called process fidelity, a fidelity based distance measure between two operations, whereas 
here we investigate in the completely bounded (CB) distance or gate trace distance, a trace- 
distance based measure. Thus, our method is expected to contribute to the problem of 



programmable quantum processors and related subjects [3l|, l32, l33| in future investigations. 
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The paper is organized as follows. Sec. [TTl gives basic formulations and main results. We 
define the error probability in realizing the quantum NOT gate based on the CB distance. 
We subsequently show that a pure input state gives the worst error probability. This enables 
us to assume, without loss of generality, that the input state is a pure state. In preparation 
for deriving the lower bound of the error probability, in Sec. IIIII we generally describe the 
maximum trace distance between the two output states from the realization and from the 
ideal quantum NOT gate. In Sec. IIVI we introduce the conservation law into the discussion. 
By minimizing the error probability over arbitrary choices of the evolution operator obeying 
the conservation law, we give a lower bound which depends only on the ancilla input state. 
In Sec. El we optimize the ancilla input state and derive a general lower bound expressed 
as a function of the size (the number of qubits) of the ancilla. Chebyshev polynomials of 
the second kind, a family of orthogonal polynomials, are used to solve this problem. To 
show the tightness of the bound, in Sec. IVIt we consider classically complete realizations, 
realizations which correctly carry out the quantum NOT operation when the input state is 
a computational basis state. We obtain the attainable lower bound for classically complete 
realizations. This result also shows that the general lower bound can be attained up to 
constant factor of the ancilla size. In the final section, we summarize our study and comment 
on the direction of future studies. 

II. FORMULATION AND MAIN RESULTS 
A. Qubits and conservation laws 

The problem to be considered is formulated as follows. The main system S is a single 
qubit described by a two dimensional Hilbert space TYs with a fixed computational basis 
{|0), |1)}. The Pauh operators Xs,Ys, and Zs on Hs are defined by Xs = |0)(1| + |1)(0|, 
Ys = -i|0)(l| +i|l)(0|, and Zs = |0)(0| - We refer to Xg as the quantum NOT gate. 

We suppose that the computational basis is represented by the z-component of spin, and 
consider the constraint on realizing the quantum NOT gate Xg under the angular momentum 
conservation law. More specifically, we assume that the control system is described as an A^- 
qubit system A also called the ancilla, and that the interaction between S and A preserves 
the z-component of the angular momentum of the composite system S + A, and study the 
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unavoidable error probability in realizing the quantum NOT operation. 

Each qubit Aj for i = 1,2, - ■ ■ ,N in the ancilla A is described by a two dimensional 
Hilbert space TYa^- Accordingly, the Hilbert space TYa of the ancilla A is the tensor product 
Ha = ^fLi'HAi, and the Hilbert space Ti. of the composite system S+A is = 7is®'^A- The 
observable Zs on Hs is identified with 2's®/ai ®/a2 ® ' ' -^/AiviWhere Ja^ for i = 1, 2, ■ ■ ■ , 
is the identity operator on TYa^, respectively. Let Z^^ be the Pauli Z operator on Ti-Ai, which 
is also identified with the corresponding operator on Ti. The sum of Pauli Z operators on 
A is denoted by 

N 

Za = '^ Za,, 

i=l 

and the corresponding sum of S + A is denoted by 

Z = Zs + Za- 

Let U be the evolution operator of S + A during the interaction between S and A to 
realize the quantum NOT gate on S. We assume that U satisfies the conservation law 

[U,Z]=0, (1) 

where [U, Z] = UZ — ZU. We shall show that the conservation law ([T]) causes unavoidable 
decoherence in realizing by U. 

To obtain the error probability, we describe the output state of S resulting from the 
evolution of S + A. Let ps and pa be states of S and A, respectively, so that the input state 
of S + A is the product state ps ® Pa- Then the output state £u,pj^{ps) of S is given by 

^c/,PA(Ps) = TrA[f/(ps®PA)f/^], (2) 

where TrA [■] is the partial trace over TYa- On the other hand, for the perfect quantum NOT 
gate, the output state £xsi.Ps) of S would be 

ExM = XspsXs^ (3) 

In the following sections we shall show that there exists an unavoidable error probability 
of the output state ([2]) in realizing the output state ([3]) under the conservation law ([1]). The 
unavoidable error probability for any unitary operator U satisfying the conservation law ([T]) 
will be evaluated to be at least 

- 1 — cos 

2V iV + 2/ 
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for the worst input state ps of S and for the best input state pa of A, and the achievabihty 
to this lower bound will be shown asymptotically. This lower bound is much tighter than 
the lower bound j^g^^^^^ anticipated from the previous investigations for other gates as to be 
shown numerically. 



B. Error probability and CB distance 

To state our results more precisely, we introduce the following definitions. Any pair 
{U, Pa) consisting of a unitary operator U on Hs ® Ha and a state pA on Ha is called a 
gate implementation or simply an implementation with ancilla A. Every implementation 



{U, Pa) determines the trace-preserving completely positive (CP) map Su^pj^ of the states 
of S by Eq. ([2]) called the gate operation determined by (f/, Pa); see Ref. [?] for trace- 
preserving CP maps in quantum information theory. An implementation (f/, Pa) is said to 
be conservative if it satisfies Eq. ([T]). We consider the problem as to how accurately we 
can make the gate operation Su,pj^ to realize the quantum NOT gate Sxs- The worst error 



probability of this realization is defined by the completely bounded distance 3J, l35|] (the 
CB distance, or the half-CB- norm-distance) between Su^pj^ and Sxs, given by 

= sup D {£u,pA ® (p) , £xs ® 'i'dn (p)) , (4) 



n,p 



where -D(-, ■) denotes the trace distance (or the half-trace- norm-distance) [7| defined by 

D{pi,p2) = ^Tr[|pi -P2I] 

for any states pi and p2 of S, idn is the identity operation on an n-level system E, and p 
runs over the density operators on S + E. Since the trace distance of the output states 
can be interpreted as the achievable upper bound on the classical trace distances, or the 
total-variation distances, between the probability distributions arising from any measure- 



ments on those states 



71], the CB distance can be interpreted as the ultimate achievable 



upper bound on those classical trace distances with further allowing measurements over the 



environment with entangled input states; see, for example, [36| for a discussion on the en- 
hancement of channel discriminations with an entanglement assistance. Thus, we interpret 
DcB {£u,pA ! ^Xs ) as the worst error probability of Su,pa realizing Sxg ■ The phrase "error 



probability" in the following discussion means the CB distance Clearly, 

DcBi£u,PA, £xs) > maxL)(£'c7,pA(ps),i^Xs(ps)), 

PS 

and minimizing maXp^ D{Su,pj^{ps)j£xsiPs)) over all the conservative implementations 
{U, Pa), we find 

-DcB {£u,pA 5 ^^s ) 

> mm maxD{£u,pAiPs),£xsips)) ■ (5) 

iU,pA) PS 

The right-hand side of this inequality can be interpreted as a precision limit of the quantum 
NOT gate under the conservation law ([1]). If the limit could take zero, it might be considered 
that there exists a perfect realization in Su,pa- However, we show that such a realization 
does not exist because of the conservation law ([1]). 

C. Sufficiency of pure input states 

Now, we shall simplify the maximization over the input state ps by showing that it suffices 
to consider only pure state ps- To show this, we use the fact that the output trace distance 
is jointly convex in its inputs: 

^ {^u,pA ( P^P^) ' ( XI P^P') ) 

i i 

< piD {£u,PA {Pi) > £xs {Pi)) , (6) 

i 

where Pi = 1 and pi > 0. This follows easily from the joint convexity of the trace distance 
and the linearity of operations Sxs and Su,pA- 

From the above inequality, a pure input state certainly gives the maximum of the trace 
distance. To see this briefiy, let ps = Xli l"^*) ("^jI' where Xli — ^ 1i — 0- Then, 
there exists a pure state such that 

i i 

<Y,q,D{8u,PAmm.^x,mm)) 

i 

<D{Su,pAm{i^3\)^£xsm{m- (7) 

Thus in considering maXpg D{£u^pj^{ps),Sxs{Ps)), we shall assume in later discussions with- 
out loss of generality that the input state ps is a pure state. 
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D. Pure conservative implementations 



An implementation (U, pa) is said to be pure if pA is a pure state. In this case, we shall 
write {U,pa) = {U, \A)) if pa = \ A){A\. In the following sections, we shall mainly consider 
the case where the ancilla state is a pure state. Here, we shall show a purification method 
that makes any general conservative implementation a pure conservative implementation, 
so that every conservative implementation with qubit ancilla has a pure conservative 
implementation with [logg rank(pA)l qubit ancilla, where rank(pA) denotes the rank of 
Pa- 

Let {U, Pa) be a conservative implementation with qubit ancilla A. Then, we have the 
spectral decomposition 

R 

Pa = 5^Pj#i)(</>il, (8) 

where R = rank(pA), {4'j\'t'k) = Sjk, Pj > for all j, k = 1, . . . , R, and J2jPj = 1- Let A' be 
the A^' qubit ancilla system extending A satisfying A^' = A^ + [log2-R]. Let \A') G Ha' be 
such that 

R 

where G Ha'^a, {^j\^k) = ^jk for all j,k = We define a unitary operator U' on 

® Ha ® Ti-A'-A by t/' = f/ ® /, where / is the identity operator on Ha'-a- 
Now, we consider the implementation [U' , \A')). It is easy to see that U' satisfies the 
conservation law [U', Z] = 0, where Z is the sum of Pauli Z operators in S + A'. We shall 
show the relation 

£u,pA = £u',\A')- (10) 

Let ps be any input state. Then, by Eq. (Q we have 

TrA'-A[ps^\A'){A'\] = ps®PA. (11) 

We also have 

£u',\A')iPs) 

= Tta'[U'{ps^\A'){A'\)U'^] 

= TtaTta'-aHU ® I){ps ® \A'){A'\){U^ ® /)] 

= TrA[UTTA'-A[ps®\A'){A'\]U^]. 
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From Eq. (fTT]) . we have 

£u',\A'){ps) = TrA[f/(ps <^ Pa)U'']. 

Since ps is arbitrary, Eq. ffTOj) follows from Eq. ([2]). 

The implementation (f/', is a pure conservative implementation and has N' = N + 
\\0g2 rank(pA)l qubit ancilla. 



E. Gate fidelity and gate trace distance 

For any two trace-preserving CP maps Si and S2 their distance measures are defined as 
follows. The gate fidelity [3| F{Si,S2) between Si and S2 is defined by 

F {Si, S2) = mfF (Slips), S2{ps)), (12) 

Ps 

where ps varies over all the states of S and F{-,-) in the right-hand-side denotes the fidelity 
defined by 

Fipi,p2) = TT[{py'p2pyY'] (13) 

for all states pi and p2 of S. By the joint concavity of the fidelity 0, P- 415] the infimum in 
Eq. (fT2|) can be replaced by the one over only all the pure states ps of S. 
We define the gate trace-distance D{Si,S2) between Si and S2 by 

D{Si,S2) = snpD{Si{ps),S2{ps)), (14) 

PS 

where ps varies over all the states of S. By the result obtained in subsection III C[ the 
supremum in Eq. (HM can be replaced by the one over only all the pure states ps of S. 
For any state pi and any pure state p2, the fidelity and the trace distance are related by 

D{pi,p2) >1-F{pi,p2f 

(see Eq. (9,111) of Ref. j3]). Since Sx^iPs) is a pure state provided that ps is pure, we have 

D{Su,pM^^xM) >l-F{Su,pM^£xMf 

(15) 
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for any pure state ps of S. Taking supremum over all the pure states ps of the both sides 
of Eq. f[T5l) . for any implementation {U, pa) we obtain 

> l-F{Su,p^,8x,r. (16) 

In Ref. 0, the realization of the Hadamard gate Hs = (l/y2)(|0)(0| + |1)(0| + |0)(1| - 
has been considered and it has been proved that for any pure conservative implemen- 
tation (f/, \A)) with N qubit ancilla A, we have 

1 

AN~+A'' 



l-F{£u,\A),£Hsr>,,^^, (17) 



where Shs{ps) = HspsH^ 39]. Since any conservative implementation (U, pa) with 
qubit ancilla A can be purified to be a pure conservative implementation {U', \A')) with 
+ [log2 rank(pA)l qubit ancilla A', we have 

Since [log2 rank(pA)l < 2A^, we conclude that every conservative implementation [U, pa) 
with qubit ancilla A satisfies 

l-FiSu,pA,£Hsr>Y^j^,- (19) 

In other words, we have 

min^,„ax|l-f £„3)^1>J^. (20) 

where [U, \A)) varies over all the pure conservative implementations with qubit ancilla 
A, and we have 

mmma.|l-F(f„,..£„.n>^^. (21) 

where (f/, pa) varies over all the conservative implementations with A^ qubit ancilla A. 



F. Main results 

Unfortunately, the method for deriving Eq. (fTTI) cannot be applied to the quantum 
NOT gate. In this paper we develop a new method for analyzing the gate trace-distance 
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D(Su,pj^,Sxs) instead of considering the gate fidelity F{£u,pj^,£xs) and we shall prove the 
following relations. In section |Vl we shall show that any pure conservative implementation 
(?7, \A)) with qubit ancilla satisfies 

1 / 2tx 



D{Su,\A),Sxs) > 2i^-"°'ivT4v'- ^^^^ 
It follows from the above, any conservative implementation {U, pa) with N qubit ancilla 
satisfies 

D(Sup.,Sxs) > -(l- cos— — — ^ -). (23) 

\ u,pA^ xsj - iV + log2rank(pA) +4/ ^ ^ 

An implementation {U, pa) is called a classically complete implementation of the quantum 
NOT gate, or classically complete implementation for short, if it satisfies 

W(|0)(0|) = |1)(1|, (24) 
= I0)(0|. (25) 

In section IVIt we shall consider classically complete pure conservative implementations. We 
shall find the attainable lower bound for this case, so that we obtain 

mill max D{Su,\A) {ps),Sxs (Ps)) 

(U,\A)) PS 

1 / 27r \ 

= ol-cos— -, (26) 



2V A^ + 2. 

where {U, \A)) varies over all the classically complete pure conservative implementations 
with N qubit ancilla A provided N is even, and we obtain 

min max D {Su,\a) (ps) , (Ps) ) 

(C/,|A» PS 

1 / 27r \ 
= - 1 - cos ) , (27) 

provided is odd. From the above, any classically complete conservative implementation 
{U, Pa) with A^ qubit ancilla satisfies 

g(g....g.v.)>^(l-cos ^^,_^^^^^;^^^^^^J . (28) 

Since A^ + log2 rank(pA) < 2A^, from the above we have 

1 r 27r 
1 — cos 



2 V A^ + 

> min maxL)(%p^(ps),^Xs(ps)) 

{U,pa) PS 

> -fl-cos^^y (29) 
- 2 V A^+ 1/' ^ ^ 
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where (?7, Pa) varies over all the classically complete implementations with qubit ancilla. 
From Eqs. ([22]) and ([27]), we have 

1/ 27r 
- 1 — cos 



2 V + 1 

- min maxD(£t/,|A>(ps),^Xs(ps)) 

(U,\A)) PS 

1 / 27r \ 
> - 1 1 - cos — ) , (30) 

- 2V iV + 4/' ^ ^ 

where {U, \A)) varies over all the pure conservative implementations. Finally, from Eqs. ( 123|) 
and fl27|) . we have 

1 / 27r 
- 1 — cos 



1 



> mm maxD{Su,pj^{ps),Sxs{ps)) 

{U,pa) PS 

> -fl-cos— (31) 
- 2V iV + 2/' ^ ^ 

where {U, pa) varies over all the conservative implementations with qubit ancilla A. 



III. LOWER BOUND OF GATE TRACE DISTANCE 

In this section, we investigate the maximum trace distance over all possible input states 
of S for given U and pa in a general way without considering the conservation law. 



A. System input state and trace distance 

We start with a description of the output states controlled by any unitary operator U on 
Hs ® Ha- Any pure input state {ip) of S can be described as 

\i;)=a\0)+P\l), (32) 

where |ap + = 1. We assume that the input state of A is a pure state \A), so that the 
input state of the composite system S + A is the product state {tp) ® \A). When |0) or |1) is 
an input state of S the corresponding output state of S + A can be generally expressed as 



um®\A)) 


= |0)g 






U{\1)®\A)) 


= |0)g 


D \Al) 4 


-\1)®\A\), 



13 



where 1^4*) e TCa for i,i = 0, 1. Normahzing these states gives 



S)ir+iii^?)f = 1, 

i)r+\\\A\)r = 1. (34) 



The output state of S + A corresponding to can then be expressed as 

U{\^)^\A)) = «(|0)®|A°) + |1)®|A?)) 

+(3 {\0) ® \Al) + \1) \A\)) . (35) 

Normahzing Eq. fl35|) gives 

Re[a*P{{A'^,\Al) + {A'M\))]=0- (36) 

The output state £^[/j^)(|?/')) := £^;7ja)(|V') (^1) of S is given by the partial trace of Eq. fl35|) 
with respect to A as follows. 



SuMm = TTA[Um{^\®\A){A\)U^] 

= (i«riii^^)ir + + a*p{A',\Ai) + mwAimio) (oi 

+ (H'iA'M'o) + c^f3*{Al\A'o) + + \P\'{A\\Al)) |0) (1| 

+ + + + \P\'{Al\Al)) |1) (0| 

+ (i«riii^?)ir + + + m\Ai)r) \i) (ii . (s?) 

On the other hand, if the quantum NOT gate were to be perfectly realized, the output state 
Sxsili^)) ■= ^Xs{\'^){'^\) would be given by 

Sx^m = Xs 1^) (^1 Xl = \(3\'\0) (0| + a*m (1| + a(3*\l) (0| + \a\'\l) (1| . (38) 

We now consider the trace distance between Su^\a){\'^)) and SxH'ip))- Note that the trace 
distance between two-dimensional states, and a'', can be described as 

Di^', ^'') = ^\4i-<\'-i4o-<oMi-0, (39) 
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where afj = {i\(T''\j) for k = ^,77. Using this relation, the trace distance 
D{Sxs{m,£u,\A){\m is 

D{£u,iAm),£xs{\m 

-0/^i'-(i«nii^o)ir+«/5*(^^i^^)+«*/5KiA^) + i/3niiA^)ir)] 

x[\a\'-{\an\A'i)r + aP*{Al\A'l) + a*(3{A'i\^^^^ (40) 



Let 60 = II 1^0) f andei = |||Ai)f . Then |||A?)f = I-60 and \\\Al)r = l-e, by Eq. dMl). 
Thus Eqs. ([36]) and (gOD give 



D{Su,\A)m,Sxs{\m = {|«*/5(1 - KI^S)) - «/5*(Al|A°) - |anA?|A°) - m'{Al\A^o)\ 

+ [- \a\'eo + |/?rei - 2Re {a* P{A',\Al))Y]\ (41) 
Clearly [{-\a\'^eo + l/Spei) - 2Re (a*/3(^°|^J))]^ > 0, and hence we obtain 

Di£u,iAm),£xsi\m 

> \a*(3 (1 - {A'ilAl)) - aP*{A\\Al) - |ap(A?|A°) 

-m'{Al\Al)\. (42) 

B. Lower bound for maximum trace distance 

In the following, we shall prove that for any U and \A), we have 

maxD {Su,\A){ps),SxsiPs)) > ^ |l - Kl^o)| , (43) 

by considering the maximization of Eq. (1421) over the input state lip) of S. This means that 
the output trace distance must satisfy Eq. fH3|) for any interaction and any input state of 
A. 

The proof is as follows. We consider the input state {ip') = a\0) + (3\1) such that |a| = 
|/5| = ^. Let e be such that = ie^^ and < < 2n. Then Eq. (02]) gives 



1 
2 



> he'' (1 - {Al\Al)) - e~''{A\\A',) - (A?|A°) 



{A\\Al). (44) 
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Here three complex numbers, 1 — {Ai\Al), —{A\\Aq), and — (A^IAq) — (AJ|Aq), which are 
determined by U and \A), can be expressed as 

l-{Al\Al) = ne'^\ 
-{Al\A',) = r2e'^^ 
-{A^K) - {Al\Al) = r,e^^\ (45) 

where n > and < (pi < 2t\ for i = 1, 2, 3. Then ri = |1 - and 

> i|^^+^2e'(-''-^^+'^^) + r3e'(-'+'^^-*^)|. (46) 

Note that Eq. fj46|) is maintained for any Q which is independent of U and \Al) . Hence, we 
consider the following two cases. In the first case, suppose that U and 1^4) satisfy r2 > rs. 
In this case, for the input state of S with d = |(02 — 0i), we have 

Di£u,\A)im).£xsim)) 

1 , 

> 2 1^1 + ^2 - 
1 

> -ri. 

- 2 

Thus, there exists a state of S that satisfies D(Sxsi\4'))i£u,\A)i\4'))) > l^^i in the case 
where r2 > r^. In the second case, suppose that U and \A) satisfy r2 < r^. In this case, for 
the input state with 6' = 03 — 0i, we have 

DiSu,\A)m),sxsm)) 

> ^\ri +^3 -r2| 
1 

> -ri. 

- 2 

Thus, there exists a state of S that satisfies D^Su^ia)^^)) , Sxsil"^))) > in the case 
where r2 < r3. We therefore conclude that for any U and \A), there exists a state \ip) of S 
such that the input state ps = satisfies 

D{Su,\A){ps),Sxs{ps)) >l\l- {AllADl (47) 
16 



This completes the proof. 

In Eq. fj43|) . if the inner product (A'JIAq) could take one by a certain choice of U and 
\A), the lower bound could take zero. This may mean a perfect realization of Su^\a) exists. 
However, we will show in the following sections that the inner product cannot take one by 
assuming the conservation law ([T]). This result will give us a precision limit of the quantum 
NOT gate. 

IV. PRECISION LIMIT GIVEN THE ANCILLA STATE 

In this section, we derive the lower bound which depends on the input state of the ancilla 
system by minimizing the right-hand-side of Eq. fHTl) over the evolution operator U under 
the conservation law. 

A. Constraints on ancilla input states 

We start with the description of the input state of A. The sum of the Pauli Z operators 
on A is the operator Za on TYa given by 

N 
i=l 

We denote the eigenspace in Za of an eigenvalue ^ by E^^. The eigenvalues are — 2n, 
where n = 0,1,2, N. The dimension of the eigenspace of the eigenvalue — 2n is 
dn = (jvi^)!n.r Note that the Hilbert space of A is the direct sum of the spaces Eff^^n fo'^ 
n = 0,l,-- - ,7V: 

Ha = ©to<-2n- (48) 

Therefore, for any input state \ A) of A there exist G C and G E^^2n with || || = 1 
satisfying 

N 
n=0 

Normalizing Eq. ( 149|) gives 

N 

I]l«n|' = l. (50) 

n=0 

17 



Next we describe the output state of S + A under the conservation law. Let E^^ be the 
eigenspace of an eigenvalue m = 1, — 1 of Zs, and be the eigenspace of an eigenvalue A 



of Z, where Z 

where = 0, 1, ■ ■ 
of the space E^^ 



Zs + Za., which has 
A 



2n, 



(51) 



A^, + 1. The eigenspace E^ can be expressed by the tensor product 



E 



Za 



i^_2n s-iid the space E^l ® E^^2n follows: 



7?^ 



e!^®eI^, 

{E^-®El^^;)®{E^J,®Elt2): 



E' 



E' 



N+l 



{E^-®Et%)®{E^.\®Et%+2) 



E^_\ ® E^^ 



(52) 



Note that the conservation law ([T]) can be equivalently expressed by the relation 40 1 

UEi C £;f (53) 

for all A. Eqs. ( 1521) and (155]) then show that the output state U{\0) ® \4>n)) an element of 
the subspace (Ef^ ® ^^^2n) ® (^-f ® ^J-2(n-i)) forn = 1, 2, ■ ■ ■ , TV, since 

um®\<l^n)) e ^(i?^®i5^^J 



C UE 



E^^ ^ 



©(£;!f ®Ej-,(„„^). 

Similarly, the output state U{\0) ® |0o)) element of the subspace -Ef® 

t/(|o)®|0^)) e f/(E^®E^-) 



(54) 



-E^^, since 



C UE 



N+l 



Za 



(55) 
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Therefore, by Eqs. iH) and there exist {{(p^fo) G Ef/^^^ and G £^^l2(n-i) 

such that 

t/(|o) ® 10;^)) = |o) ® 1(0;^)°) + |i) ® (56) 

where = 0. Normahzing Eq. ( l56l) gives 

iii(0n)^)ir+iii(0ii)?)ir = i- (57) 

Similarly, for the output state ?7(|1) ® there exist |(0;^+i)o) G -E^l2(n+i) 

and 1(0;^)!) G 

-^^-2n s^*^^ that 

u{\i) ® 10;!)) = |o) ® |(0;!^i)^) + |i) ® |(0^)i), (58) 

where |(0;^+i)o) = 0. Normalizing Eq. fl58l) gives 

lll(0n+l)o)ir+lll(0n)})ir = l- (59) 

We can now obtain useful relations for the output state of S + A under the conservation 
law. For the output state f/(|0) O |v4)), Eqs. 1^ and (IMD give 

N 



um®\A)) = \o)0(j2<^n\{<i>X) 



n=0 
N 



(E^-K^^-i)?))- (60) 



n=0 



Similarly, for the output state U{\1) |A)), Eqs. (09]) and ([58]) give 



N 



Ui\l)®\A)) = \0)^(Y,an\{<Pti)l) 



n=0 
N 



+ (E^-K'^n)!))- (61) 



n=0 



Comparing Eq. fl33|) with Eqs. fl60|) and fl6T]) . we obtain the following relations: 



AT 



\<) = E«"i(<^n)o)' 

n=0 

Af 

n=0 

Af 

n=0 
N 

= J2^n\{<P^)\). (62) 



n=0 
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B. Optimization of gate trace distance by ancilla input 



We can now estimate the inner product {Al\Al). By Eq. ( !62|) . 



N 



(A'Ml) = E <^n*aM^,_,mti)l)^ (63) 



n,n'=0 
lA \0\fAA \1\ 



where the inner product ((0„/_i)i|(0„+i)o) is given as 

((0n'-l)?l(0ll)o) 

for n' — 1 ^ n + 1, 

{i<PtiM<Pti)h) for n'-l = n + l. 



Therefore, 



N~2 



n=Q 

l((0?l(Oo)l < lll(0?)lllll(Oo)ll<l- 



Thus, 



N-2 



(64) 



Ki^o) = (65) 

n=0 

By the triangle inequahty, we have 

N-2 

\{Al\Al)\ < E K^,\ K\ mti)i\i<Pti)l)\- (66) 

n=0 

From Eqs. fl5U]l . flFTI) . and flM?]l . the Schwarz inequahty gives the relations 

N-2 

|an+2||an| < 1, (67) 



(68) 



< E|a„+2||aJ < 1, (69) 

n=0 

SO that the maximum of KA^I^q)! is at most '^^=0 I 1 1 c^n I • Therefore, the minimum of 
i|l — (A^I^q)! in the right-hand side of Eq. fH3l) is at least i(l — J2n=o Wn+2\Wn\)- Since 
in the above argument the unitary operator U was arbitrary but satisfied the conservation 
law, we have 

minmaxZ:' {Su,\a){Ps),Sxs{Ps)) 

U PS 

N-2 



n=0 



20 



where U varies over all the unitary operators on 'Hs^'Ha satisfying Eq. ([T]). This is a useful 
inequality that allows us to evaluate a lower bound of the quantum NOT gate given the 
input state of the ancilla system. For example, if constant, such as 

for all n = 0, 1, ■ ■ ■ , A^, then whatever evolution operator is used, an error probability 
determined by Eq. ( ITOi) is unavoidable. 

The following questions regarding Eq. fl70p still remain: What is the lower bound over 
the input states of the ancilla system? Can we reduce the lower bound to zero by choosing 
appropriate input states of A? In the next section, we will quantitatively show that there 
exists a non-zero lower bound of the error probability for any input state of the ancilla system 
and any evolution operator. In order to obtain the bound, it is necessary to minimize Eq. (170|) 
over the input states of A under condition ( 1501) . 



V. PRECISION LIMIT GIVEN THE ANCILLA SIZE 

We consider the maximization of J2n=o kn+2||on| over input states of the ancilla system 
to minimize the right-hand side of Eq. fl70|l under condition fl50|) . In the first place, we 
show that this problem can be reduced to the derivation of the maximum eigenvalue of a 
symmetric matrix. Secondly, we explain how to derive the maximum eigenvalue, making use 
of the recurrence formula of Chebyshev polynomials of the second kind. We finally describe 
the lower bound of the quantum NOT gate which depends only on the size of the ancilla 
system. 



A. Lower bound and eigenvalue problem 

The summation J2n=o kn+2||«n| can be divided into two parts, the summation of odd 
subscripts, such as |ao| |a2|, |a2| |a4|, ■ ■ ■ , and that of even subscripts, such as |ai||a3|, |a3||a5|, 
• • ■ . For even A^, 

Af-2 



n=0 



JV-4 JV-2 
2 2 



|fl2r+l| |fl2r+3| + | «2r 1 1 fl2r+2 1 , (72) 



r=0 r=0 
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where > 2. For odd N, 

N-2 

|'2„+2| jOnl 

n=0 

JV-3 JV-3 
2 2 

= |Q2r+l| |Q2r+3| + | Q2r 1 1 Q2r+2 1 ; (73) 

r=0 r=0 

where > 3. We now assume that is even for simphcity; we will comment on the case 
of odd later. To rewrite the summation, we define an (A^ + l)-dimensional vector by 

A"^ = [ |ai|, lasi, . . . , |aAr_i|, |ao|, |a2|, • • • , \aN\], (74) 

where the odd indexed (resp. even indexed) elements are in the first (resp. second) half 
elements of the vector, and the number of those elements is y (resp. y + 1). The summation 
can then be expressed by a matrix and the vector A as 

N-2 
n=0 





1 



■■■ 

1 


■•• 





At 




1 



■■• 





■■■ 





1 
1 






■■• 







1 




where the matrix has four submatrices. The upper left (resp. lower right) submatrix is the 
Y X Y (resp. (y + 1)x(y + 1)) matrix with all the first sub diagonal entries one and all 
the other entries zero. The upper right (resp. lower left) submatrix is the y x ( y + 1) (resp. 
(y + 1) X y ) matrix with all the entries zero. Taking the complex conjugate of both sides 
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of Eq. 



gives 



N-2 



Wn+2\\ 



n=0 



At 





1 
1 



■■ 

to Eq. 



Therefore, adding Eq. 

Af-2 



n=0 



At 





1 









6]) gives 








1 



2 





2 


1 

2 





1 



i 



A. 



I 



i 



A, 



(76) 



(77) 



where the upper left and the lower right submatrices are symmetric with all the first sub- 
diagonal and superdiagonal entries 1/2 and all the other entries 0. Let At^^ and At^^^^ be 
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two vectors defined by 



"-odd 













l«2|, 


a4 , . . 


• ) \(^n\ ] , 



(78) 



and S"; be an / X / symmetric matrix defined by 



i 

1 i 



S, 



i 

■•. 



1 

■ 2 

I 



Then, Eq. fl77|l can be written as 

N-2 
n=0 



< IIAoddr -Si^ + IIAc 



S N 



(79) 



(80) 



where si is the maximum eigenvalue of the symmetric matrix Si. Recall that HAoddP + 
IIApvpiiP = A"!" ■ A = 1, and thus 



max 



'N~2 

|an+2| |On| 



n=0 



max 



(81) 



where the maximization in the right-hand side means selecting the larger of Siv and Si 



-1- 



Taking the difference between Eqs. (1721) and (1751) into account, we apply the same analysis 
for odd A^. Then, we have 



max 

E|an|2 = l 



'N-2 



n=0 



S N+l . 

2 



(82) 



In this way, the maximization of the summation Yln=o kn+2| kn.| under condition ( l50l) re- 



duces to the derivation of the maximum eigenvalue of the symmetric matrices Sn_ and Si 



+1- 



B. Eigenvalue problem and orthogonal polynomials 



Next we shall determine the maximum eigenvalue, as mentioned above, and give the lower 
bound of the quantum NOT gate. It is well-known that the eigenvalues and the eigenvectors 
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of the matrix Si are obtained from a recurrence formula of orthogonal polynomials as follows 



StI . |38 |. Chebyshev polynomials Wi{x) ioi I = 1,2, ... oi the second kind are defined by the 

sin(/ + 1)9 



relation 



Wi{cosi 



sin 6 



(83) 



where < 6 < n, and are polynomials of the precise degree I, and satisfy the recurrence 
formula 



xWoix) = ^W,{x), 

xWiix) = ^Wi+i{x) + ^Wi^i{x), 

where / > 1. The roots x = xi^k of the equation Wi{x) = is given by 

kiT 



xik = cos 



/ + 1 



for k = 1,2, . . . ,1. Let W'''(x;^fc) be an /-dimensional vector defined as 

Wt(x;,fc) = [Wo{xi,k), Wi{xi,k), ■■■ , Wi.,{xi, 
Since Wi{xi^k) = 0, Eqs. ([85]) and give 



SiWix, 



i 

1 i 

i ••. 







1 

■ 2 

1 



Woixi^k] 

WAxLk) 



Wi-i{xi^k) 



lWi{xi,k) 
IWoixi^k) + \W2{xi,k) 

^Wj.iixi^k) + ^Wj+i{xij 



Xl,k'^{xi^k)- 



(84) 
(85) 

(86) 

(87) 



(88) 



Thus, the vector W(x;^fc) is an eigenvector of Si with eigenvalue x^k- Therefore, the maxi 
mum eigenvalue of Si is 

TT 



Si = Xil = COS 



l + l 



(89) 
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and the corresponding eigenvector is given by 



' . (j + l)7r 

sm — ; — 

l + l 



TX 



sm 



/ + 1 



j=0 



(90) 



C. Derivation of lower bound given the size of ancilla 



We have found the maximum eigenvalue, and thus we can now describe a lower bound of 
the error probability in realizing the quantum NOT gate. For even A^, Eqs. ( IHTi) and ( 189|) 
give 



N-2 

max } |a„+2| knl = cos 
1 ^-^ 



J2M 



n=0 



27r 

iV + 4' 



(91) 



Recall that the minimization of Eq. (ffOl) over the input states of A is derived from the 



maximization of J2n=o I '^■"+2 1 l^-nl- Thus 



mm maxD{Su,\A){ps),Sxs{ps)) 

(U,\A)) PS 



2tx 



cos 



(92) 



Similarly, for odd 



mill max L)(6:[/,|A>(ps),^Xs(ps)) 

{V\A)) PS 

1 / 2n 

> - 1 — cos — 

- 2V iV + 3 



(93) 



Here cos is greater than cos-^^^, and hence we have finally obtained the lower bound 
for the error probability of any realization of the quantum NOT gate with A^-qubit control 
system under the angular momentum conservation law as 



mill maxD{Su,\A){ps),^Xsips)) 

{U,\A)) PS 

1 / 2n 

> - 1 — cos 

- 2V iV + 4 



(94) 



for any A^(> 2). The bound depends only on the size of the ancilla system: the larger N, 
the closer to zero is the lower bound. 



According to previous works 12|, ll9j based on the uncertainty principle, it may be ex- 
pected that the lower bound of the quantum NOT gate scales with the inverse of N as 
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FIG. 1: Plot of the lower bounds as a function of N. The solid line shows the lower bound 
^(1 — cos j^^) of the quantum NOT gate in Eq. (f94l) . The dashed line shows the lower bound 



4(jv^+i) Pi'sviously obtained for the Hadamard gate in Ref. [19l |. 

4(Ar2-i-i) ~ However, the new bound has the leading order |(1 — cos -j^^) ~ so that 
the lower bound obtained here is really tighter than that as depicted by Figure [H 



D. Lower bound: general case 

We have considered the case where the ancilla state is a pure state. In the following 
we shall consider the general case. Let (f/, pa) be a conservative implementation with 
qubit ancilla A. Then, its purification {W, \A')) is a conservative pure implementation with 

+ [log2 rankpA] qubit ancilla A' such that Su^pj^ = Su',\a')- Applying Eq. (!94l) to Sui^\a'), 
we have 

max D {£u,p^ [ps ) , ^Xg (ps ) ) ) 

Ps 

>-fl-cos— -V (95) 

- 2V iV + log2rankpA + 4/ ^ ^ 

and from + log2 rankpA < 2A^, and we conclude 

min max D {Su.p^ (ps ) , ^Xs (ps)) 

{U,pa) PS 

> ifl - COS— (96) 
- 2V iV + 2/' ^ ^ 

where {U, pa) varies over all the conservative implementations with qubit ancilla. 



27 



VI. LOWER BOUNDS FOR CLASSICALLY COMPLETE IMPLEMENTATIONS 
AND THEIR ATTAINABILITY 

In the preceding section, we have shown that a general lower bound for the error 
probability in realizing the quantum NOT gate is given by the 1 — cos(l/A^) scale for 
the ancilla size A^, instead of l/N"^ scaling already known for some other gates. Since 

2[1 — cos(l/A^)] = l/N"^ — 1/(12A^^) H , the new scale has the same leading order as 

up to constant, but it is natural to ask if the higher order terms are really meaningful. Here, 
we shall answer this question, so that the 1 — cos(l/A^) scale is the best result. To show 
this, we shall show the attainability of a lower bound with the 1 — cos(l/A^) scale for classi- 
cally complete conservative pure implementations. Thus, a classically complete conservative 
implementation exists even with only 2 qubit ancilla, whereas the substantial error occurs 
when the input state is a superposition of computational basis states. This result also shows 
that the general lower bound for conservative implementations with qubit ancilla can be 
reached by a classically complete conservative pure implementations with 2N qubit ancilla. 

A. Classically complete pure implementations 

Let (f/', \A')) be a classically complete conservative pure implementation. Then, we have 
the following relations 

U'i\0)®\A')) = 

U'{\1)(^\A')) = \0)®\A'l), (97) 

where \A'1) and \A'l) G Ha- 

First, we discuss the constraint on the input state \A') of A imposed by the above 
relations. To illustrate this, we describe \A') as 

N 

\A')=Y.a'J<P^'), (98) 

n=0 

where \(f)n') are normalized vectors in the eigenspaces E^^2n all n = 0, 1, ■ ■ ■ ,N, and we 
have ^^=0 I'^'nl'^ ~ ^- Suppose that the input state of S is |0). Recalling that relation fl55|) 
holds by the conservation law, the output state corresponding to the input state |0) ® |0o') 
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can be written as 

f/'(|0) ® 10^')) = e^<^'|O)®|0^'), (99) 

where e'^'^' is a phase factor. Thus the output state corresponding to the input state |0) ® \A') 
can be expressed as 

f/'(|0) ® \A')) 

N 

= a'oe^^'|O)®|0^') + E<f^'(|O)®l<^"'))- (100) 



n=l 



Comparing with Eq. flW|) . must be zero. Similarly, a'^ must be zero, considering the input 
state 

We now describe the output state in S from {U', \ A')) for any pure input state \ip). This 
is given by the partial trace of the output state in S + A with respect to A: 

= TrA [U'iltp) ® \A')) m ® (A'l) U'^' 

= \Pm{0\ + a*P{A",\A'l)\0){l\ + aP*{A'l\A",)\l){0\ 

(101) 

Here, we use abbreviation such as S^ip)) ■= S{\ip){ip\) for any operation S. The trace 
distance between the ideal quantum NOT operation fl55]) and £^;7/ja')(I^)) is then 

DiSxsm),su',\A'){\m = ii - (102) 

Thus, the derivation of the lower bound for the gate implementation {U', \A')) can be reduced 
to estimating the maximum value of {A'1\A'l), which is very similar to the general analysis 
of Sec. IIVI However, this case differs from the general analysis in that = oat = 0. Taking 
this condition into account, \A'1) and \A'l) can be written as 

N-l 



n=l 
N-l 

\<) = E<l('^"+i)o)' (103) 



n=l 



where \{4>n'-i)i) \{(f)n'+i)o) normahzed vectors in the eigenspaces -E'^^2(n-i) ^^"^ 
-^^^2(n+i)' respectively. Thus, 

N-3 

\{A'l\A'l)\ < (104) 



n=l 
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and therefore, 

min max D {£u',\a') (Ps ) , £xs (Ps ) ) 

U' Ps 

> 2(l-El<+2lKl)- (105) 

n=l 

Since the discussion in Sec. |V]can be apphed to minimizing Eq. (11051) over the input states 
of A, we see that for even 

{U',\A')) PS 'I /\r / 

1 / 27r \ , , 

> - 1-cos— ). 106 

- 2V N + 2J ^ ' 

This lower bound is shghtly larger than the one for the general case; the difference comes 

close to zero for large N of the ancilla system. We shall comment on the odd case later. 



B. Attainability of the lower bound for classically complete pure implementations 

Next we show that there exists a classically complete implementation ([/', \A')) which 
attains the lower bound |(1 — cos;^^^). We begin by describing the input state \A) as 
follows. Let |(e„)*) be fixed orthonormal bases in eigenspace Eff^^n 

\{enf)Men)').--- Menf)r-- MenY-). (107) 
for n = 0, 1, ■ ■ ■ , A^, where dn = ^[(N-ny. • 1^ addition, AJ^^^ and Al^^^ are two vectors: 

^Idd = [oi, «3, Os, ■ ■ ■ , aAT-i ] , 

Alven = [ 02,04,06, ■■■ ,aAr-2]. (108) 

where A^^^ (resp. Aj^^j^) is a ^ (resp. y — 1) dimensional vector whose entries are indexed 
by odd (resp. even) numbers. We assume that these vectors satisfy 

Aodd = 7^W(xiv J, 

2 

Aeven = 0, (109) 

where Cn = [W{xN^yyV{xNj^)]^. It follows that ||AoddP = 1 by normalization. We 
assume that the input state \A) is given by 

|i) = ^Sj(e.)^). (110) 

n=l 
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Recall that W(xjv ^) is an eigenvector with the maximum eigenvalue of Sn_. Then the 
coefficients a„ satisfy the following equation: 

7V-3 
n=l 

= Sn 

2 

= cos-^. (Ill) 

Constructing the evolution operator U can be accomplished by determining the transfor- 
mation for all orthonormal bases. We require that U satisfy the following conditions. For 
n = 1, 2, ■ ■ ■ , AT, 



U{\0)®\{en)')) = |1)® |(e„_i)^), 
f/(|l)® |(e„_i)i)) = |0)® |(e„)i), (112) 

and for all bases except those that appear in Eq. flll2p . 

t>(|O)®|(e„)0) = |0)®|(ej^), 

l(en)^)) = |(en)T (113) 

These requirements determine one-to-one mapping on the orthonormal basis, {|0) ® 
I (e„)*) , 1 1) ® I (e„,)*) }, and hence there uniquely exists a unitary operator U fulfilling the above 
requirements. Note also that U satisfies the conservation law ([1]), since from Eqs. (11121) and 
(11131) we have the relations UE^ C for all A, which are equivalent to the conservation 
law, as seen in Eq. (1531) . 

We now describe the output state of (f/, \A)) and the trace distance between the ideal 
output state and that of {U, \A)). The output states for |0) and |1) can be generally written 
as 

i/(|o)® |i)) = |o)® |i°) + |i)® 

U{\1) ® \A)) = |0) ® \Al) + |1) ® \A\), (114) 
respectively, where \A^,) G Ha with i,j = 0, 1. On the other hand, by the definitions of U 
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and \A), we have 

f/(|0)® |i)) = t/(|0)® ^a„|(eO^) 

n=l 

N-1 

n=l 



n=l 

7V-1 



|0)® (^a„|(e„+i)i)). (115) 



n=l 



Thus we have the following relations: 



1^^) = 0, 



iV-l 



n=l 
Af-1 



n=l 



|i}) = 0. (116) 

Let i^(7j^)(|V^)) be the output state of S from (f/, \ A)). The trace distance between Sxsili^)) 
and ^e expressed in the same way as for Eq. fj4T|) so that we have 

D{£a,iA)im,£xs{\m 

= {\a*p[l-{A',\Al))+aP*{A\\A',) 
-\a\'{Al\A',)-\P\'{A\\Al)\' 
+ [{-\a\% + \P\%) 

-2Re{a*(3{A',\Al))Y}\ (117) 

where |||^o)|P = cq, = ii. However, in this case, eo = ei = from Eq. (11161) . and 

therefore 

DiSo,\A)im,SxM))) = |a*/?(l-(i;|ij))|. 
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FIG. 2: Distribution of |a„| with odd subscripts for N = 100 which gives the lower bound in 
Eq. (jl06p . This figure shows -^r—Wn-i (xn ^) as a function of odd n. 



' N 



Recall that Ke^)^) are orthonormal bases. Then, Eq. (11 111) gives 

AT-l 

= anan'((e„_i)^|(e„'+i)^) 

n,n'=l 
N-3 



n'=l 



271 



COS ■ 



Thus, 



Di£^^^^^{m,£xM))) = |a*/3(l-cos^ 



2tt 



Since the right-hand side is maximized where 



^, we have 



(118) 



maxD{S^,^M'^)),Sxsi\'4'))) 

W 

1 / 27r \ , , 

= - 1-cos— . 119 

2V N + 2J ^ ' 

That is, the model (t/, 1^4)) attains the lower bound in Eq. (11061) . Notice that our model 
(JJ , \A)) has a distribution of |a„|, as given by Eq. (I109p . Figure [2] describes the distribution 
for = 100. From a qualitative point of view, in order to reduce the lower bound of the 
quantum NOT gate, an input state of the ancilla system should be prepared which has a 
sufficiently thick distribution in the neighborhood of eigenvalue 0, rather than a constant 
distribution, such as that given by Eq. ( 1711) . 

For odd A^, the lower bound can be given by setting the input state and the evolution 
operator as those analogous to the case of even A^. The bound is |(1 — cos;^|^). The 
attainability of this bound is also proved by the analogous argument. 
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Thus, we have shown that 



min max D {Su,\a) (ps ) , ^Xg (ps ) ) 

(f/,|A>) PS 

1 / 27r \ , , 

o 1-cos— — , (120) 



2 V N + 2 

if N is even and 

min maxD{Su,\A) (ps) , ^Xs (Ps) ) 

{U,\A)) PS 

if is odd, where {U, \A)) varies over all the classical complete pure implementation with 
qubit ancilla. 

For arbitrary A^, we conclude as a common lower bound 

min maxD{Su,\A){ps),£xs{Ps)) 

{U,\A)) PS ^ 'I /V / 

> - 1 - cos— , 122 

- 2V iV + 2/' ^ ^ 

where ([/, \A)) varies over all the classical complete pure implementation with N qubit 
ancilla. 

We have considered the case where the ancilla state is a pure state. The lower bound 
for the general case is obtained by the previously developed purification argument, and we 
conclude the following relations. We have 

max D {£u,pj^ (ps ) , £^Xs (Ps ) ) 

PS 

1 2n 



> - 1-cos ; , (123) 

- 2V A^ + log2rankpA + 2/ ^ ^ 

for any classically complete implementation {U, pa) with N qubit ancilla, and 

min max D {Su,p^ (ps) , ^Xs (ps ) ) 

(.U,Pa) Ps 

>U^-^os-!^), (124) 



2 V N+l 

where {U, pa) varies over all the classically complete implementation with qubit ancilla. 
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VII. CONCLUDING REMARKS 



In this paper, we have studied the precision hmit of the quantum NOT gate or the bit flip 
gate, one of the most basic gates in quantum computation, represented on the single-spin 
computational qubit by considering the angular momentum conservation law obeyed by the 
interaction between the computational qubit and the control system supposed to comprise 
many qubits. Actually, we have considered the effect of the angular momentum conservation 
law only in the direction same as the computational basis, usually set as the z direction. 
Then, the conserved quantity and the computational basis are represented by the Pauli Z 
operator, whereas the quantum NOT gate is represented by the Pauli X operator. Thus, it 



is expected that this non-commutativity leads to a precision 
In the previous method which was used for other gates 



12 



imit of the gate operation. 



19| , one finds a way in which 



the gate under consideration is used as a component of a measuring apparatus, applies the 
quantitative generalization of the Wigner-Araki-Yanase (WAY) theorem to this measuring 
apparatus, and obtains the lower bound of error probability. For the Hadamard gate, one 
finds that it is used to convert the Z measurement to the X measurement, and that Z 
measurement can be done without error under the conservation law of the z component. 
Then, one can conclude that the inevitable error of the X measurement, calculated from the 
quantitative version of the WAY theorem, is yielded from the converter using the Hadamard 
gate. This and similar arguments cannot be applied to the quantum NOT gate, since the 
quantum NOT gate does not convert the direction of measurement, but simply flips the 
measured bit. 

In this paper, we have developed a new method for obtaining the inevitable error probabil- 
ity by evaluating the maximum trace distance between the output from the gate realization 
and the output from the ideal gate. The previous method naturally leads to a lower bound 
for the infidelity (one minus the squared fidelity). Since the infidelity is dominated by the 
trace distance, the new method gives a tighter lower bound for the error probability. 

The new method is based on a straightforward evaluation of the trace distance of two 
output states, and enables us to find the precision limit Eq. (170]) . explicitly described by 
the input state of the ancilla system. It is thus possible to obtain information on how much 
an ancilla input has an inherent error probability in itself. The correspondence between the 
two methods is not easy to elicit, but it is an interesting problem for future studies that 
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would lead to a deeper understanding of precision limits to quantum control systems. 

We have also obtained the lower bound flMl) expressed by the size of the ancilla system, by 
minimizing Eq. (1701) over the input states of A, using Chebyshev polynomials of the second 
kind. The lower bound is much tighter than the scaling expected from the previous result 
based on the WAY theorem. Since the quantitative generalization of the WAY theorem has a 
close relation to the universal uncertainty principle for measurement and disturbance 19|, |20| , 



the previous lower bound for pure conservative implementations is based on the variance of 
the ancilla state, and scales as ~ whereas the new method revealed the lower 

^ as a tighter bound. The higer order terms in |(1 — cos;^^^) 



bound |(1 



cos 



2n ^ 



is considered to be meaningful, since the lower bound |(1 — cos j^^) is attained among 
classically complete pure conservative implementations. Interestingly, the attainability result 
shows that the best ancilla states to attain the lower bound are not maximum variance states, 
nor uniformly distributed states, but those states with the distribution determined by the 
recurrence relation characterized by Chebyshev polynomials. 

Although our study has assumed that the ancilla system consists of N qubits for com- 
parison with the previous research, the present method is not restricted to this particular 
control system, and it can be readily applied to other control systems, such as atom-field 
systems, where the present method would lead to a lower bound that scales as the inverse 



of the photon number 



of programmable quantum processors [28|, l29l 
investigations. 



231 ]. Our method will be also expected to contri 

and related subjects 



Dute to the problem 
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331] in future 
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